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We consider a thick brane model supported by two interacting scalar fields in 7D and 8D general 
relativity. Using the special type of a potential energy, we obtain numerically the regular asymptot- 
ically flat vacuum solutions. A possibility of obtaining the similar solutions for an arbitrary number 
of the extra spatial dimensions is estimated. 



I. INTRODUCTION 

The idea about the multidimensional Universe appeared in the 1920s for unification of gravitational and elec- 
tromagnetic interactions. Fast progress of the quantum theory has overshadowed researches of multidimensional 
gravitational theories unifying gravitation and other types of interactions. Just in the 1980s a great deal of interest 
has been resumed within the framework of superstring theory pretending to unify all of the fundamental interactions 
in one theory. Perhaps the most important problem in the multidimensional theories is non-observability of extra 
dimensions. Firstly, it was proposed that extra dimensions are compactified and unobservable up to the Planck energy 
scales. Despite the numerous attempts for solving the problem of compactification of extra dimensions, one can think 
that the satisfactory solution of the problem was not found at the moment. 

At the end of the 1990s, the idea that extra dimensions could be noncompactified (infinite), but the matter is 
somehow concentrated on a 4D surface (brane), was proposed [1|] (see also the earlier works in this direction Q). On 
the basis of this idea the new ways for solution of some old problems of high-energy physics (the problem of mass 
hierarchy, stability of proton, etc.) were proposed. At the moment, the models of 4D-branes embedded in 5D and 6D 
spaces Q and also crossing 5D-branes in 6D Q are under consideration. 

It is supposed in brane models that various types of matter (bosons, fermions, gauge fields and so on) are localized 
on a hypersurface (brane) embedded in an external multidimensional spacetime. However, in contrast to the multi- 
dimensional Kaluza-Klein theory, extra dimensions could be macroscopic ones and, generally speaking, non-compact 
ones. Another essential difference consists in different behavior of gravitational and matter fields: if the former ex- 
ists and propagates in the bulk, the matter fields are localized on a brane and they are 4D objects. At such an 
approach, however, the situation could be realized when a multidimensional gravitational field could be also localized 
on a brane 0, @] (see also Q), i-e. effectively it becomes four dimensional despite the extra dimensions could be 
macroscopic ones. The arising effective 4D gravitational constant is defined not only by physics on a brane but by 
influence of extra dimensions also. These effects could be directly checked in experiments and observations both on 
small and cosmological scales (see, e.g., Ref. Q). 

A special interest consists in consideration of various scalar fields within the framework of the scenarios of the world 
on a brane. It is caused by that scalar fields are widely used in particle physics, astrophysics and cosmology H| . Within 
the framework of brane theory, scalar fields were considered in different aspects in Refs. [1, HO, HH 112, El, 13, 13 (fo^ 
a review, see jl6|). This interest to the scalar fields is explained also by relative simplicity of equations and solutions 
obtaining in models with using of them. It allows to make qualitative analysis of equations and find sufficiently clear 
physical interpretation of results by analogy from other fields of physics. 

One of the main problem at consideration of brane models is that in almost all brane models one uses infinitely 
thin branes that, of course, is absolutely unsatisfactory from the physical point of view. The point is that obtaining of 
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self-consistent solutions of gravitational equations interacting with matter and describing a thick brane is not a simple 
mathematical problem. At the moment there exist the following thick-brane solutions: in Ref. [l3| such solutions are 
obtained as monopole-like solutions. The set of gravitating scalar fields with non-trivial topological configuration far 
from the brane is considered. It was shown that such solutions exist for co-dimension of the brane more or equal to 
2. Unfortunately, for co-dimension equal to 1 there is no topologically non-trivial configurations and, consequently, 
there is no thick-brane solutions in this context. Also in Ref the monopole-like solutions with some potential of a 
scalar field with broken symmetry were considered. The regular or periodic solutions for a global monopole with an 
arbitrary number of the extra spatial dimensions n and for a local monopole with n = 3 were obtained. In Ref. [l9j 
a thick-brane solution in non-Riemannian generalization of the 5D Kaluza-Klein theory was obtained. In such an 
approach, affine connection, except Christoffel symbols, has also the Weyl term which is described by a scalar field. 



In Refs. [20|, [2l| the thick-branes were obtained at special choice of a potential of a scalar field. It has allowed to 



integrate once the field equations for the scalar field. The essential factor of obtaining of these solutions consists in 
existence of the potential of the scalar field unbounded below, i.e. V{(j)) —oo as </) ^ ±oo. 

In the light of aforesaid, it becomes interesting to consider various types of multidimensional brane models with 
scalar fields for the purpose of finding out of possibility for existence of localized solutions with finite energy for 
different number of extra dimensions. In this paper we show that two nonlinear gravitating scalar fields could create 
a 4D brane in 7D and 8D spacetimes. From the physical point of view the situation is the following: an interaction 
potential of these fields is constructed in such a way that it has two local and two global minima. It means that there 
are two different vacuums. The multidimensional space is filled by these scalar fields which arc located in that vacuum 
in which the scalar fields are in the local minimum. There is a defect in the form of a 4D brane on the background 
of this vacuum. 

For 5D case such a model was considered in Ref. [l^], for 6D one - in Ref. [1^. The obtained solutions indicate a 
possibility of localization of scalar fields of special type on a brane. 

The search of thick brane solutions is very complicated mathematical problem. Indeed, only three type of solutions 
are known till now: (a) thick branes with a special choice of a potential of a scalar field [20|; (b) thick branes with 
topologically non-trivial scalar fields 17] and (c) thick brane with two nonlinear scalar fields [23] ■ Thus the goal of 
our paper is to extend the solutions found in Ref 's to higher dimensions. We do not consider trapping of any 

matter on such branes because it is an independent and complicated problem demanding a separate consideration. 
We hope to investigate trapping of matter in our future investigations. 



II. GENERAL EQUATIONS 



In general case of 13 = 4 + rt dimensional gravity the action can be written as follows \24 



S = / d^X\/D 



-R + L„ 



(1) 



where M is the fundamental mass scale and n is a number of extra dimensions. As a source of matter fields we 
chose two interacting scalar fields <p, x with the Lagrangian: 



where the potential energy V{(p, x) is: 



V{^, X) = ^(^' - rnlf + ^{X'- mlf + ^'x' V,. 



(2) 



(3) 



(This potential energy was obtained in [25j at approximate modeling of a condensate of gauge field in the SU(3) 
Yang-Mills theory.) Here and further capital Latin indices run over A, i? = 0, 1,2,3, and small Greek indices 
= 0,1,2,3 refer to four dimensions; Ai,A2 are the self-coupling constants, mi, 1112 are the masses of the scalar 
fields ip Xi respectively; Vq is an arbitrary normalization constant which will be chosen below from physical reasons. 

Use of two fields ensures presence of two global minima of the potential ^ at (p — 0,x — ± TO2 and two local 
minima at x = 0, = ± mi for values of the parameters Ai,A2 used in the paper. The conditions for existence 
of the local minima are: Ai > 0,mf > A2TO2/21 and for the global minima: A2 > 0, > Aim^/2. The presence of 
these minima has allowed to find solutions localized on the brane for 5D and 6D cases in Refs. [l^], [1^ when the 
solutions have tended asymptotically to one of the local minima. 

Variation of the action ([T]) with respect to the D-dimcnsional metric tensor qab led to Einstein's equations: 
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where i?^ and are the I?-dimensional Ricci and the energy-momentum tensors respectively. The corresponding 
scalar field equations could be obtained from ^ by its variation with respect to the field variables ip, x- Then these 
equations will be: 
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dV 



(5) 



III. 7D CASE 



In this case n = 3 and the metric could be chosen in the form [24| : 

ds^ = (P^{r)T]apdx"dx'^ - X{r)(dr^ + r^d^l), 



(6) 



where r]ap is a flat 4D spacetime metric and the metric functions and A depend only on the extra coordinate r. 
dl72 = d9^ + sin'^ 6dip^ is an angular part of the metric depending only on 6-th and 7-th coordinates. 
Then, using (H)),© and (O, one can get the system of gravitational equations in the following form: 
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where a prime denotes differentiation with respect to r and the corresponding components of the energy-momentum 
tensor are: 
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Multiplying the equation ([9]) by 3/4 and subtracting its from ([7]), and also subtracting the equation ([8]) from 
have: 
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where we have used the following rescalings: r r/M^^^, ip — > AI^^^ip, x ~^ M^^^x^ ^1,2 
equations are supplemented with two equations for the scalar fields which are obtained from ([5 
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2X1^^' = ^^ [2x' + Ai(¥''-m2)] , 
X" + (7 + 4^ + Ij) X' - Ax [2p' + A^ix' - ml)] 



mi 2 . These 
with use of ([6]): 

(14) 
(15) 



The set of equations (|T2)) - (|T5)) is a set of nonlinear equations with solutions whose behavior depends essentially from 
values of the parameters mi, m2, Ai, A2. As it was shown in Refs. [l^l, [111, specifying some values of the self-couphng 
constants Ai, A2, a problem of search of regular solutions for systems similar to (IT^ - (IT5]) reduced to evaluation of 
eigenvalues of the parameters mi, m2. Only for these values of the parameters the regular solutions with finite energy 
exist. 

We will search for such regular solutions in a whole space. As it is not possible to find analytical solutions, we 
search for a solution numerically. In order to make a solution regular everywhere, let us consider a solution near the 
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brane, i.e. at r « 0. Let us search for this solution as a series: 
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where the subscript denotes a value of the function at r = and the subscript 2 denotes a second derivative of the 
corresponding function at r = 0. Substitution of these series into the equations (fT2|l - (fT5l) gives the following values 
for 02, A2,(;32,X2: 

A2 = ^A^y(O), 02 = -^0oAonO), (20) 
15 15 

V2 ^ ^Xofo [2x1 + M {vl ~ mj)] , X2 = ^AoXo [2¥'o+A2(xo-"i2)] • (21) 

The obtained equations (|12p - p3|) could not be solved analytically, and so we solve them numerically. But also at 
a numerical analysis there are some difficulties. The point is that the preliminary numerical analysis shows that a 
solution regular everywhere exists not for all mi 2 but for some special values of this parameters. This means that we 
have deal with a nonlinear eigenvalue problem. The functions f{r),x{''') ^-re eigenf unctions, and the parameters 77ii_2 
are eigenvalues, and the condition Ai 7^ A2 should be satisfied. In particular, we chose Ai = 0.1, A2 = 1.0. 



A. Numerical solution 



In this section we describe a method of numerical solution of the equations (|T2)) - p5l) in details. As it was mentioned 
above, the special feature of the set of equations is that determining of the functions tp and x is a nonlinear eigenvalue 
problem. Usually, numerical solution of one ordinary differential equation for search of eigenfunctions and eigenvalues 
is carrying out by the shooting method. The point of the method is that one, using a method of step-by-step 
approximation, try to find with necessary accuracy some eigenvalue at which an eigenfunction is regular. For example, 
in such a way one can find a discrete energy spectrum of a particle in a ID potential well of an arbitrary shape. 
Unfortunately, this method does not work in a case when one has a eigenvalues and eigenfunctions problem with two 
variables. It takes place in our case: we have two eigenfunctions ip,x E^nd two eigenvalues mi^2- 

Therefore we will search for a numerical solution in the following way. Searching for eigenvalues of the equations 
P^ - (fT5)l . we solve on each step only whether the equation (fH|) or Since there is also another function in this 

equation, its value takes from a previous step. Having found (with some accuracy) some values of the functions 
(f, X ^^^d eigenvalues mi. 2, one should insert them into the Einstein equations (|12p - (jl3p for determining of the metric 
functions 0, A. Then these functions should be inserted in one of the field equations (fT4|) - p^ which is solving as 
a eigenfunctions problem for the functions (p,x and eigenvalues toi_2 once again. These iterations repeat necessary 
number of times for obtaining of a solution in the range r S [A, rj], where r/ is some final point on the axes r. 

Note the following features of the above method of a solution search of the differential equations (fT^ - P^ : 



• One could not start a numerical solution from r = because there are terms like y{r)/r in the equations 
(|12p - (jl5p . where y{r) is one of the functions 0, (/3,x,A. Therefore one should start a numerical solution from 
r — A and use the following boundary conditions: 
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• One should not try to find a numerical solution for large r/ because any numerical method has some inaccuracy. 
This inaccuracy defines (by some unknown way) the interval [A,r^] in which an exact solution differs from 
numerical one insignificantly. 

• An eigenvalues problem is very sensible to calculated eigenvalues. For example, let us consider the equation 

fj^^iyixMxf-e]. (26) 

This equation, being considered as an eigenvalues problem with the eigcnfunction y{x) and eigenvalue e, has 
an exact solution y — tanh(a::),e = 1. If one solves this equation numerically with initial conditions which are 
followed from the exact solution 

y(0) = 0, y'(0) = l, (27) 
then one can see that for x « 6 the numerical solution differs significantly from the exact one. 
Thus we chose the following algorithm of numerical solution of the equations (fT2|) -(fT5 l) : 

1. A choice of a zeroth approximation for the function ^p{r) or x(r). 

2. Numerical solution of the equation (|14p or by the shooting method with evaluation of the functions tp{r) 
or x(r) and mi or m2. 

3. Numerical solution of the equation ([T5| or (fT^ by the shooting method with evaluation of the functions x(r) 
or ip{r) and m2 or mi. 

4. Reiteration of the steps 2 and 3 necessary number of times. 

5. Substitution of the functions tf{r) and x(r) into the Einstein equations IT^ and and their solution by a 
usual numerical method. 

6. Substitution of the functions 0(r) and A(r) into the equations or (fTS)) and reiteration of the steps 2-5 
necessary number of times. 

7. Check of the obtained solutions by solution of the initial equations (fT^ - (fT5)) with use of the calculated on the 
steps 2-3 parameters toi^2- 

For these calculations we use the program package Mathematica for numerical solution of the differential equations. 



B. Results 



Using the procedure described in the previous section, let us find a self-consistent numerical solution of the system 
([T^ - ([T5)) . We chose the following boundary conditions at r = 0: 

ipo = 73, = 0, 

Xo = VfTe, x'q = 0, 

00 = 1.0, cb'o = 0, (28) 

Ao = 1.0, = 0, 

and Ai = 0.1, A2 = 1.0. The arbitrary constant Vq is chosen as Vq = (A2/4)m2 for zeroing of the energy density as 
r — !■ 00. 

In this case regular solutions exist at mi w 2.31505626 and m2 ~ 3.08288116. The obtained solutions are presented 
in Figs. [U [5) The corresponding plot for the energy density is shown in Fig. [31 As one can see from the last figure, 
the energy density tends asymptotically to zero as r —^ 00, i.e. the scalar fields are trapped on the 4D-brane. 

It is difficult to see an asymptotic behavior of the metric functions from Fig. [2l For clarification of this question, 
let us find asymptotic solutions of the equations (fT4|) - (fT5l) in the form: 

(p^mi-S(p, X = Sx, (29) 
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FIG. 1: The scalar fields 'p{r),x{''') for the boundary conditions 
given in (|28(l . The solid lines correspond to the 7D case, the 
dashed lines - 8D case. 



FIG. 2: The metric functions (j>{r),X{r) for the boundary con- 
ditions given in (I28|l . The solid lines correspond to the 7D 
case, the dashed lines - 8D case. 




FIG. 3: The energy density Tq (r) from (|10[) . The solid line corresponds to the 7D case, the dashed line - 8D case. 



where Sip,Sx ^ 1 as r — > cx3. Such the asymptotic behavior of the fields corresponds to tending of the solutions to 
the local minimum aX cp = mi,x = 0. Then the right hand sides of the equations P^ - (fT5|) tend to zero, i.e. one 
has the following asymptotic equations for the metric functions: 

A" 3/A'y 12 4 / 0' 13A'\ 12 0'A' 



Let us find solutions of these equations in the form: 



= 0. (31) 



(32) 
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(The subscript "oo" indicates the asymptotic value of the variable, Ci > 0, C2 > and a,/3 are some arbitrary 
constants.) One can check that in the case under consideration a and (3 are equal each other. Then we have from 

« 0OO - — , (34) 
r 

X ^ X^ + ^. (35) 
r 

Inserting these solutions in (|14p - (fl31) and taking into account we have: 

Sip" + ^S(p' = 2X^AimjSip, (36) 
Sx" + l5x' = X^{2ml-A2ml)Sx (37) 

exp (-A/2AooAim^ r) 

Sif « ^ ^, (38) 

r 

exp (-VAoo(2toi - A2TO2) ^) 

<5x ~ C;, ^ ^, (39) 

r 

where C,^,Cx are integration constants. Thus one can see from ([29]) . ([32|) . (p3|) . ([38|) and p9|) that the asymptotic 
behavior corresponds to a 7D Minkowski spacetime with the scalar fields energy density equal to zero. 



with regular solutions 



IV. 8D CASE 



Now let us consider the 8D problem. For this case n — A and we chose the metric as follows: 



<j)^{r)rji,^dx''dx'' - A(r) {dr^ + [dO^ 



(40) 



here /i, = 0, 1, 2, 3 refer to four dimensions; r, 9, ip, ^ are extra coordinates; 0, ?/;, ^ are polar angles on a 3D sphere; 
77^1, — {+1, —1, —1, —1} is the 4D Minkowski metric. Inserting the metric (|40p in (H])-®, one can obtain the following 
equations: 
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(41) 
(42) 
(43) 
(44) 
(45) 



We will search for solutions of this system by analogy with the 7D case. Using the same boundary conditions 
one can obtain the results presented in Figs. [T][3l 

It is interesting to estimate an asymptotic behavior of the solutions as r ^ 00. For this purpose let us search for 
a solution of the equations (|44| - ((45|) in the form (|29|) . Then the right hand sides of the equations (j4ip - (|43|) tend to 
zero and we will search for solutions for the metric functions in the form One can obtain from (|4ip - (P51) 

that a = (3 = 2, i.e. we have the following asymptotic behavior of the metric functions: 

_ (46) 
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Inserting these solutions in (|44p - ([45l) and taking into account ([29l) . we will have: 



with regular solutions 



+ = 2AooAim?(5(^, (48) 
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5^'+'U^ = XW^ml-k^mDSx (49) 
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where C<^, C;^ are integration constants. Thus one can see from (|29p . (j46|) . (|47p. ([50|l and ([STjl that, such as in the 7D 
case, the asymptotic behavior corresponds to a 8D Minkowski spacetime with the zero scalar fields energy density. 



ARBITRARY NUMBER OF THE EXTRA SPATIAL DIMENSIONS 



Obtaining of the numerical regular solutions for the 7D and 8D cases allows us to hope that similar solutions could 
exist and for an arbitrary number of the n extra spatial dimensions. Unfortunately, one cannot obtain numerical 
solutions for the arbitrary n because it is not possible to eliminate n from the equations. However, we can estimate 
a possibility of obtaining of such solutions. Let us use for this the Einstein equations obtained in Ref. [l^l for the 
generalized Z?-dimensional metric 

ds^ = (j)^{r)'n^p{x'')dx°'dx'^ - X{r){dr^ + r'^dnl_i) , (52) 

where dV,'^_i is the solid angle for the (n — 1) sphere. It is convenient to rewrite the Einstein equations from [24] in 
the form: 
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As one can see from the equations ([55 )) - ((55|) . the left hand sides could be regular at r = if the boundary conditions 
A'(0) = OjJ>'{0) = are hold true. Just these conditions for the metric functions were used by us at consideration of 
the 5D [22] and 6D |23| cases and also the 7D and 8D problems in this paper. Apparently, the situation remains the 
same and for a case of an arbitrary number of the n extra spatial dimensions. 

On the other hand, consideration of an asymptotic behavior of the metric functions (f>, A for the 7D and 8D cases 
shows that their behavior could be described for the arbitrary n > 6 in the following form: 

« -^oo - (58) 
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Such a behavior corresponds to fast transition of the solutions to a Minkowski spacetime. At the same time, an 
asymptotic behavior of the scalar fields could be described as follows: 



exp 



^ ~ mi-C^ , (60) 



exp 



X - ^(^rzi)7^ • (61) 

So we have the regular solutions for the arbitrary n both in neighborhood of zero and as r — > cx3. Note that it is 
possible to expand the solutions near zero in a Taylor series (see ([T6|) - (fT9| ) and obtain solutions with any necessary 
accuracy. One might expect that these solutions could be joined smoothly and they would be regular in the entire 
range < r < oo. 



VI. CONCLUSIONS 



We have shown that two gravitating nonlinearly interacting scalar fields could form a 4-dimensional thick brane 
configuration in 7D and 8D spacetimes. Consideration of the problem reduced to investigation of the system of 
ordinary differential equations. It was shown that the regular solutions with the finite energy density exist only for 
some values of the masses of the scalar fields mi,m2 at some self-coupling constants Ai, A2 given beforehand. I.e. 
the problem of a search of eigenvalues of the parameters mi,TO2 for the system of nonlinear equations ([4])- ([5]) was 
solved. There are regular solutions with an asymptotically flat metric (l){oo) = (/)oo, A(oo) = X^o and with the zero 
energy density of the matter fields at mi k, 2.31505626 and TO2 ~ 3.08288116 (for the 7D case) and mi « 2.25005 and 
1712 ~ 3.115 (for the 8D case). In both cases the solutions start with the boundary conditions ([^5]) with subsequent 
transition to the local minimum of the potential ([3]). 

A numerical analysis shows that the solutions exist since the potential V{ip, x) has the local minimum. Our attempts 
to find a solution tending to the global minimum as r ^ 00 were unsuccessful. The possible reason is the Derrick's 
theorem [26j which forbids an existence of static regular solutions with finite energy for scalar fields in dimensions 
more then 2. At numerical solution (on the step 4) we obtain a solution which is regular both in a whole 7D and 8D 
spacetimes. At first sight this solution seems to be forbidden by the Derrick's theorem but more careful analysis shows 
that the obtained solution avoids the conditions of the given theorem because the scalar fields tend asymptotically to 
the local minimum but not to the global one. 

Let us discuss a question about stability of the obtained solutions. If one consider these fields as classical ones, then 
it is possible to test the stability by the following way: (a) first of all, one can investigate small perturbations by a 
standard method; (b) then one should investigate large perturbations; (c) following which it is necessary to consider 
a behavior of this solution after quantizing of the scalar fields. In the last case the quantum fields could tunnel from 
a region of the local minimum (false vacuum) to a region of the global minimum (true vacuum). But in this case, 
apparently, the solution should decay into some waves because the existence of a static solution with scalar fields 
tending asymptotically to a global minimum is forbidden by the Derrick's theorem. 

But situation could turn out to be more interesting: the point is that in Ref. [25| there are some arguments in 
favour that the scalar fields considered in this paper are quantum non-perturbative condensate of SU(3) gauge field. 
Briefiy these arguments consist in the following: components of SU(3) gauge field could be divided by some natural 
way into two parts. The first group contains those components which belong to a subgroup SU{2) e 51/(3). The 
remaining components belong to a factor space SU{3)/SU{2). Non-perturbative quantization is carrying out by the 
following approximate way: it is supposed that 2-point Green functions could be expressed via the scalar fields. The 
first field If describes 2-point Green functions for SU(2) components of a gauge potential, and the second scalar field 
X describes 2-point Green functions for SU{3)/SU{2) components of a gauge potential. It is supposed further that 
4-point Green functions could be obtained as some bilinear combination of 2-point Green functions. Consequently, 
the Lagrangian of SU(3) gauge field takes the form 

At such an interpretation, a question about stability of the obtained thick brane solution becomes rather nontrivial 
problem. In this case the obtained solution describes some defect in a spacetime filled by a condensate of a gauge 
field. A question about stability demands a non-perturbative quantized consideration of SU(3) gauge field with Green 
functions depending on time. Difficulty consists that the suggested in Ref. (25| approximate method of description of 
Green functions could be used only in a stationary case. 

Finally, in section |V] we have discussed a possibility of obtaining of thick brane solutions with an arbitrary number 
of the n extra spatial dimensions. It was shown that there are regular solutions for regions near zero and as r ^ 00 
for the arbitrary n. It allows to hope for existence of solutions joining these regions smoothly. 
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The whole set of resuhs obtained in works '22*1 for 5D and [l^ for 6D cases, and also in this paper, allows us 
to speak about possibility in principle of localization of the scalar fields with the potential ([3]) on a brane in any 
dimensions. One has every reason to suppose that existence of similar regular solutions is possible for models with 
greater number of extra dimensions. In future we suppose to investigate models with arbitrary number of extra 
dimensions for clarification of this question. 
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